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1. Introduction 

Ever since the seminal work of Broadbent and Hammersley in 1957 jl], the lasting 
interest in directed percolation and related statistical models has been nurtured by 
several motivations. Firstly, very anisotropic fractal structures are observed in Nature 
in a variety of situations where a bias is introduced by an external field or a gradient, 
ranging from river basins [g] to stress cracks in metals |§ and to lightning discharges. 
The preferred direction may alternatively correspond to time, with fractal patterns 
appearing in the dynamics of the systems considered. Pomeau thus proposed that 
in simple models of turbulence the spatio-temporal structure of the active regions could 
be described by directed percolation (DP in the following), which was later confirmed 
by extensive numerical simulations || [7J. More generally the question arises if the 
large-scale properties of the observed structures depend on the details of the dynamical 
processes occurring on short scales, or if they are more universal and can be described 
in terms of purely statistical models. 

Another motivation is that directed systems may be studied using renormalization 
group (RG) tools, providing a possible approach to nonequilibrium phenomena. The 
introduction of a preferred direction is a relevant perturbation in the RG sense, it is 
expected to change the universality class of a phase transition and modify the critical 
exponents in a non-trivial way. In particular the DP threshold is in the same universality 
class as Reggeon field theory ||. In spite of much work, however, exact results are known 
only for simpler variants and DP itself remains a challenge for theorists, being one of 
the few fundamental problems to remain unsolved in two dimensions - see M for a 
recent review of this model and its position among models of non-equilibrium phase 
transitions. As pointed out by Grassberger fl(| it is surprising in view of the central 



role played by DP that no detailed experimental tests of the predictions for its critical 
behaviour have been carried out. A serious difficulty is that in practice many effects 
can obscure the expected behaviour at attainable scales, so it is of interest to explore 
the crossover regime between DP and other universality classes. 

Directed lattice animals provide another class of fractal structures that may be 
relevant for oriented systems, and it has been suggested for example that they might 



describe the large-scale geometry of river networks [ 11] . The basic animal model consists 
in enumerating all connected clusters of N sites on a lattice but various extensions may 
be considered by introducing local interactions between the cluster sites. One can in 
this way obtain different statistical ensembles ranging from loopless structures (lattice 
trees) to fully compact ones. In contrast with the situation for DP several exact results 



are known for directed animals on various lattices [12|, [13], [L4], 15, |T^] - for recent results 



and a review of what is known rigorously see |T7J and references therein. 

A remarkable connection between the two types of systems has been found by 
Dhar fll8| , who showed that a model of interacting directed animals (IDA) can be 



mapped onto a directed percolation problem if certain relations hold between their 
respective parameters. For sufficiently strong attractive interactions the typical animals 
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correspond to percolation clusters in the dense phase above the percolation threshold. 
As a consequence the collapse transition of the IDA model is expected to be of the 
same nature as the DP threshold, in contrast to the isotropic case where the percolation 
threshold has been found to correspond to a higher-order multi-critical point in the 
phase diagram of interacting animals [|19|, [2(J . 

We consider in the present work the statistical properties of fully directed lattice 
animals with nearest-neighbour interactions on the square lattice, using a transfer matrix 
approach to calculate their correlation lengths on infinite strips of finite width. The 
results are analysed using a phenomenological renormalization group method adapted 
to directed systems |15| , |2lj , which allows an accurate determination of the critical 
fugacity. 

The paper is organized as follows: in the next section the IDA model and the 
main quantities of interest are introduced and we present the correspondences with 
a percolation problem and with the generating functions of various simple restricted 
ensembles. In section |3] we apply the phenomenological renormalization (PR) group 
approach to non-interacting animals and check its accuracy by comparing the results 
obtained using strips of widths up to n = 18 with the exact values or the best available 
ones. Section f| is devoted to the study of the behaviour at the collapse transition 
temperature, we determine the corresponding critical geometrical and thermal exponents 
and compare them with DP exponents. The low-temperature regime, which corresponds 
to the percolating phase of DP, is considered in section ^. In this regime the transition 
occurring at the critical fugacity is first-order and the grand partition function displays 
an essential singularity, but the finite-size PR analysis remains applicable and yields 
results in excellent agreement with Dhar's prediction for the critical line |18[]. In section ^| 
we obtain very precise values for the critical fugacities of various ensembles, such as bond 
animals, site or bond trees and strongly embeddable animals [ 22fl , and we compare them 
with published values when these exist. In section |7| particular attention is paid to 
the crossover region in the immediate vicinity of the collapse transition. Though the 
singularity is very weak we are able to obtain a direct estimate of the crossover exponent. 
We present analytic arguments which show that this crossover exponent is equal to the 
Fisher exponent a for DP ||23|| , in agreement with our numerical results. 



2. The interacting animal model 

2.1. Directed animals 

A directed animal is a connected cluster of N sites, some of which belong to the root, 
such that any other site of the animal can be reached from the root by a walk which never 
goes opposite to the directed axis (see figure |]). In the simplest case, non-interacting 
animals, one associates the same weight x N to each cluster having N sites, so the basic 
generating function is 

S (x) = J2 n N x N , (1) 

N 
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Figure 1. A directed site animal on the square lattice, containing N = 33 sites and 
P = 37 pairs of nearest-neighbour sites (filled circles), on a strip of width n — 4 with 
periodic boundary conditions. The root sites are denoted by an extra circle. 



where is the number of directed connected clusters having a single site (the origin) 
as root. 

The problem of enumerating such objects on lattices was recognized as an 
interesting combinatorial problem and received much attention over the last two decades, 
in part because it turned out to be possible to obtain many exact results concerning the 
fijv, whereas for undirected animals only bounds are known. In particular the generating 
function ([!]) is known for the square lattice studied here JL4], [L7|] : 



S (x) 



1 // 1 + x , 
2^1 - 3x' 



1/2 



(2) 



For physical applications the metric properties of directed animals are perhaps yet 
more important. They are best described in terms of two independent characteristic 
lengths with asymptotic behaviour of the form 



R A 



c 2 N v 



(3) 



where R" provides a measure of the average distance between the root and the most 
distant sites along the directed axis, while R L corresponds to the average distance 
between two extreme animal sites in the perpendicular (transverse) direction. The 
value of the perpendicular exponent is known to be v 1 - = 1/2 in d = 2, from the exact 
solution for various lattices [|14], pl| and in agreement with a field theory argument 
connecting the large-scale behavior of directed animals and the dynamics at the Yang- 
Lee edge singularity [T3, |T3|. The longitudinal exponent v» has proven much more 
difficult to obtain and is known only numerically |24], |25 . 



2.2. Branched polymers and interacting animals 



Interacting animals are often used as simple lattice models to describe the behaviour of 
branched polymers in dilute solutions. Indeed, field theory calculations |26] indicate that 
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the existence of loops in animals is irrelevant in the RG sense, so one expects them to 
belong to the same universality class as branched polymers and to share the same large- 
scale behaviour. A simple way to take into account the interactions with the solvent 
is to associate an energy e with every pair of neighbouring occupied sites. This energy 
corresponds to the difference between the monomer-solvent and the monomer- monomer 
interaction energies (one can distinguish between chemically bonded pairs and others in 



mere physical contact |p7| , but the present simple model already contains much of the 
physics). At an absolute temperature T the Boltzmann thermal weight associated with 
each pair is then 

w = exp(e/k B T), (4) 

where k B is the Boltzmann constant. 

A partition function for polymers made of exactly N monomers may be defined: 

Z N {T)=^{N,P)w P , (5) 
p 

where Q(N, P) is the number of clusters of N sites having P pairs of occupied nearest- 
neighbors (note that all possible polymer topologies contribute to Zn, this corresponds 
to a model of "living polymers" where the monomers may rearrange). At infinite 
temperature w = 1 and one recovers the statistics of an ensemble of equally weighted 
animals. A good solvent corresponds to e ~ (w ~ 1), a poor solvent corresponds to 
e >> (large w). The free energy per site in the thermodynamic limit is given by 

/(T)= lim - tzLlog Z N (T). (6) 

N^oo iv 

It is often convenient to consider the generating function 

G(x, w)=J2 z n(w) x n = J2 V(N, P) x N w p , (7) 

N NP 

which is analogous to a grand partition function and has a singularity for x = x c (w), 
the critical fugacity. From (§), © and (|7|) one has 

f(T) = k B T log x c (w). (8) 
2. 3. Collapse transition and percolation threshold 

When the temperature is lowered compact configurations become more probable and 
the polymers are expected to shrink from an extended state to a dense globular state. 
In the limit N — > oo a collapse transition occurs between these two states at a well- 
defined temperature (for a statistically typical polymer). This transition is described 
by a critical point, usually referred to as the #-point, which is analogous to a tricritical 
point for magnetic systems. It is believed that the average linear size of the polymers at 
the #-point scales as N UB , where vq differs from the exponents v and v c corresponding 
respectively to the extended and dense states. Various lattice models have been 
proposed for the collapse of isotropic branched polymers |27j. No complete solution 



of these models on regular lattices is known, and most of them have been studied by 
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numerical methods (see, e.g., P8[ ) such as Monte Carlo simulation, exact enumeration on 
lattices, or transfer-matrix calculations coupled with phenomenological renormalization 
(PR) @ 0. 

For directed animals the corresponding transition was first studied by Monte Carlo 
simulations [[3TJ, within an isotropic scaling approach. Although this assumption did 
not properly take into account the intrinsically anisotropic nature of directed animals, 
the value obtained for the 9 temperature is found to be in agreement with more accurate 
approaches. 

A central result, due to Dhar and rederived independently in a recent study of 



the adsorption of interacting animals [32], is the existence of a correspondence between 



the partition function of directed lattice animals with interactions between first and 
second neighbours and a site-bond percolation problem on the same directed square 
lattice. In particular the IDA first-neighbour model studied here corresponds to a simple 
bond percolation problem: a given animal A occurs with probability 

V(A) = l±l {J??-)* (I±V (9) 

and using (HD one obtains the relation 

Pr(N) = £ V{A) = I±* Z N (w) (10) 

\A\=N P 1 + " 

where the sum bears over all clusters of A" sites and 

p= l-q = (w-2)/(w- 1) (11) 

is the associated bond percolation probability. Pr(N) is the probability that for this 
value of p a percolation cluster starting from the origin contains precisely A" sites. Hence 
in the thermodynamic limit the free energy per site is given for w > 2 by 

j^f = M^-) " fon (1/A0 log Pr(iV), (12) 
Kb± 1 + q Jv^oo 



For isotropic percolation in d dimensions it has been proved []33| that above the 
percolation threshold the decay of Pr(N) is slower than exponential and that there 
exists a value p\ > p c such that for p > p± 

log Pr (AO ~ -CN l - l/d if Af^oo, (13) 

but in fact this behaviour is expected to hold for all p > p c . Physically this just means 
that most sites of a large finite cluster lie far from its external boundary and belong 
to its "bulk", whose effective volume is of the order of the geometric one, while the 
external boundary that needs to be disconnected from the infinite cluster is not fractal 
and contains of the order of N^ d ~ l '' d sites [[33]. 



In the directed case large finite percolation clusters are quite rare above p c and 
difficult to study directly, but a proof has been outlined recently that in this case too 
the decay of Pr(N) is slower than exponential |}2 and it is natural to expect that 
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log Pr(N) also follows a power law decay for all p > p c . The second term on the right- 
hand side of (0) is then negligible in the thermodynamic limit and one obtains a simple 



equation for the critical line in the low-temperature phase [18 



Xc ( w ) = ML= w 2 if w>w (0). (14) 

V 1 1+q w{w- l) 2 W V ' 

A remarkable consequence is that the free energy per site is not singular on the low-T 
side of the transition. The excellent agreement of this relation with our numerical results 
presented below confirms its validity. As a further consequence the critical temperature 
of the IDA model studied here is related to the critical probability p c for the directed 
bond percolation process on the same lattice: 

w(6) = (2- Pc )/(l-p c ). (15) 

The threshold value p c as well as the values of DP critical exponents have recently been 
estimated very accurately using low-density series expansions pi , 35, ^6 1, so the thermal 



critical exponents of collapsing directed animals are also known with high precision. 

Note, however, that the IDA model is more general than DP as it remains well 
defined for repulsive interactions (w < 1), whereas the equivalence with DP is only 
valid for w > 2. Also, the whole (x,w) plane may be given a physical meaning if the 
region x > x c (w) is interpreted as a "weak gel" stabilized through application of an 



external pressure 29 



2-4- Special ensembles 

Let us remark at this stage that the model defined by (^) contains in its phase diagram 
several special points related to ensembles defined by simple geometrical properties. 

Note firstly that for directed animals the number £ of independent loops is just 
equal to the number N 2 of sites connected to two occupied predecessors. Every such 
site contributes two pairs of occupied neighbours in (|7]) while the other occupied sites 
contribute only one, except for the origin (for simplicity we consider animals rooted at 
a single site), so the total number of pairs is P = 2N 2 + (N — 1 — N 2 ). The number of 
loops is then 

£ = P-N + l (16) 
and the loop generating function |17| is given by 

Gi(y,v) = Y,tt(N,£)y N v c (17) 

NC 

= vG(y/v,v). (18) 
In particular one obtains the generating function G st (y) of site trees (loopless animals): 
G st (y)=G l (y,0)= lim wG(x,w) (19) 

with y = wx fixed. 

A simple correspondence also exists between non-interacting directed bond animals 
(i.e. clusters of connected directed bonds) and site animals: For every occupied site 
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having only one predecessor consider the bond linking them, while for the N 2 sites with 
two predecessors consider one of three possibilities, either chosing one of the connecting 
bonds or keeping both of them. Since N2 is equal to the number £ of loops this 
correspondence implies that the generating function of these bond animals is given 
by 

G b (y) = £ Q(N, C) y N -'- c (2y + y 2 f (20) 

NC 

= ±G l {y,2 + y) (21) 

y 

= -G(x,j^-), (22) 
x 1 — x 

with x = y/(2 + y) and where we have used ([17]) and fll8|) to obtain (|2T|) and (p2j) 
respectively. 

If one now discards from the previous bond configurations all the loop-forming 
ones (i.e. those with two connecting bonds going to the same site) one is left with 
directed bond trees. The correspondence with site animals just discussed shows that 
their generating function is given by 

G bt (y) = -G l (y,2) (23) 

y 

= -G(x,2), (24) 
x 

with x = y/2. 

Finally, strongly embeddable animals |E| correspond to maximally connected bond 
configurations, in the sense that for a given site animal one systematically keeps the two 
bonds going to every site with two predecessors (bond trees correspond to minimally 
connected configurations). Their generating function is then simply 

G emb (y) = E n(N, Q y N '^ c (y 2 f (25) 

NC 

= -Gi(y,y) (26) 

y 

= G{l,y). (27) 
Detailed numerical results for these various ensembles are presented in section []. 

3. Numerical methods and phenomenological renormalization analysis 

3.1. Transfer matrix method 

Here we shall use the method of transfer matrices coupled to a phenomenological 
renormalization analysis |37 , 38| , that was applied to lattice animals in fl~5| , p9| . It is 



convenient to define the generating function Gor^ w) for the total number u;or(A^, P) 
of directed animals of sites with P pairs of nearest-neighbors connecting the points 
and R 

GW(x, w) = E ^o R (iV, P) x N w p . (28) 



NP 
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On a strip of width n with, say, periodic boundary conditions, linear recursion relations 
may be written between the u; 0R (iV, P) for two successive values of the distance R along 
the strip. The asymptotic behaviour of Gon(x, w) is dominated by the largest eigenvalue 
X n (x,w) of the corresponding transfer matrix T n : 

G 0R (x,w) ~ [X n (x,w)] R = exp ( — R - ), R^oo, (29) 

V &(x,w)J 

which defines the longitudinal (or parallel) correlation length: 

£l(x,w) = -l/\og\ n . (30) 

This correlation length diverges along the critical line x = x n (w), where x n is the 
smallest positive root of 

X n (x n (w),w) = 1 (31) 

and the generating function is infinite for x > x n (w). Let us recall that this line 
determines the free energy per site of a very large animal of width n, f n = k^T \o<gx n {w). 
In a similar way \ n (x, w) determines the average site density p n {w) of an animal on the 
strip [|T^, [2D|] through 



1 d\og[\ n (x,w)] 

Pn{w) = - 



X — X 12. 



(32) 



n d log x 

where the derivative is calculated at x n {w). 

The gel phase, for x > x n (w), will not be studied in the present work, let us just 
note that expression (|32| ) remains valid and gives the gel density |29| , but the expression 



of the correlation length is different (see for a lucid discussion of that point). 



3.2. Phenomenological renormalization in the high - temperature limit 

In the high-temperature limit (w — 1) one recovers non-interacting lattice animals 
which as mentioned above have been studied by a variety of methods (see [15, 25, 2^ 



and references cited therein). Let us recall here just two recent results. Using a 
connection of this problem to the continuous-time dynamics of a one-dimensional lattice 
gas [03 , Dhar was able to obtain a very accurate value of the longitudinal size exponent 



v" = 0.81733(5), where (5) indicates the estimated uncertainty for the last digit. This 
critical exponent was estimated more recently |E5J by the method of series expansions, 
which gave z/" = 0.81722(5). Both these results exclude the conjecture z/" = 9/11 which 
was based on phenomenological renormalization on strips of width up to n = 12 fTsU , 
and it is now generally believed that z/" is not a rational number. 

We shall first extend the latter calculation on strips of width up to n = 18. As 
we shall show, our estimate of z/" is in excellent agreement with the results just quoted, 
which provides a useful check of the overall reliability of the present approach to the 
case of interacting animals. 
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According to renormalization-group theory it is expected that for n 3> 1 and 

x — x c <C 1 the longitudinal correlation length satisfies a finite-size scaling relation 
of the form 

Ci-^Fin^ix-Xc)], (33) 

where the scaling exponent k — v^/ u 1 - and the scaling function F(z) is finite and regular 
for z = 0. This allows to obtain an estimate of the critical fugacity x c as the fixed point 
of phenomenological renormalization equations involving three different strips. Using 
three consecutive widths (n — 1, n, n + 1) the equations for the fixed point read 

£n-l{ x c,n) _ £n{ X c,n) _ £,n+l( X c,n) (34) 

(n — l) Kn n Kn (n + l) Kn ' 
where x c ^ n and K n are respectively estimates of the critical fugacity and the scaling 
exponent. Eliminating K n in (|34]) gives an equation for x Ci „: 



log = -2-2- log for x = x c , n , (35) 

a fog— d-i 



which would yield the exact value of x c if the scaling relation (|33D held rigorously. 
Finite-size estimates K n , and v\ of the critical exponents are then obtained from 

fog(d+i/ei) 



log[(n + l)/n] 



(36) 



1 _ , \og[(dd+i/dx)/(d£l/dx)} 
K log[(n + l)/n] 

v\ = n n v^ , (38) 

where all ratios and derivatives are evaluated at x = x Cjfl . 

To get accurate numerical estimates of the required derivatives we have used here 
a perturbative expansion method for the transfer matrices as described in |[40|| . This 
enabled us to compute the largest eigenvalue derivatives with an accuracy typically 
better than 1CT 12 . We also took advantage of the fact that the transfer matrix is 
singular for w = 1 (i.e., det(T n ) = 0) to reduce its size: for instance, for n = 15 the 
matrix size could be reduced from (1223 x 1223) to (198 x 198). 

3.3. Analysis of the numerical results 

In order to extrapolate these finite-size estimates two different extrapolation procedures 
were used. The first one (to be referred to as extrapolation 1) is based on a simple 
sequential fit of the form 

v£ =v ± + An~ Wl + Bn-" 2 , (39) 

where the estimates of and of A, B,lji,lj 2 can be obtained from the values of 
for five consecutive widths. One could expect that for large n these estimates take 
certain limiting values, providing the form (|39|) is well chosen. It turns out, however, 
that in some extrapolation procedures considered in this work the effective values of 
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Table 1. Critical properties of non- interacting directed lattice animals on the strips 
of a square lattice with periodic boundary conditions. The finite-size estimates 
are obtained from the renormalization group equations (|35-38). Extrapolation 1 is 



based on the three-term fit (39), while extrapolation 2 relies on the BST algorithm. 



Uncertainties in the last quoted digits are due to the extrapolation procedures and are 



shown in parentheses. 


n 


< 


v\ 




3,4,5 


0.483452612 


0.829035839 


0.333659392915 


4,5,6 


0.489970354 


0.825744198 


0.333462353775 


5,6,7 


0.493175785 


0.823887630 


0.333395848362 


6,7,8 


0.495010984 


0.822687170 


0.333367807660 


7,8,9 


0.496169426 


0.821844121 


0.333354128213 


8,9,10 


0.496951869 


0.821218766 


0.333346732444 


9,10,11 


0.497507400 


0.820736423 


0.333342415719 


10,11,12 


0.497917305 


0.820353299 


0.333339741933 


11,12,13 


0.498229180 


0.820041906 


0.333338005223 


12,13,14 


0.498472489 


0.819784080 


0.333336832461 


13,14,15 


0.498666307 


0.819567315 


0.333336014429 


14,15,16 


0.498823451 


0.819382709 


0.333335427953 


15,16,17 


0.498952800 


0.819223752 


0.333334997468 


16,17,18 


0.499060679 


0.819085572 


0.333334674963 


Extrapolation 1 


0.49999(3) 


0.81731(5) 


0.333 333 33(1) 


Extrapolation 2 


0.50000(3) 


0.81730(5) 


0.333333 34(1) 


Exact 


1/2 




1/3 



B and u>2 are rather unstable, which probably indicates that a more complicated form 
should be used for the next to leading correction. In such cases we prefer to keep only 
the first two terms on the right hand side of (^), which always seem to lead to a 
reliable extrapolation. In either case the new sequences of estimates that we calculate 
are monotonic and they depend on n more weakly than the original ones. The final 
estimates are obtained by plotting the extrapolated data versus n~ u were uj is adjusted 
to give a smooth convergence. 

The second procedure used to extrapolate for n — > oo, which will be referred to here 
as extrapolation 2, is known as the Bulirsch-Stoer (BST) algorithm [^TJ. This method 
includes a free parameter u which has to be adjusted self-consistently to minimize the 
fluctuations in the sequences of estimates produced by the BST algorithm. Let us 
note here that the ranges of uj estimated by this method are in agreement with the 
corresponding estimates for the leading correction to scaling exponent U\ from fl39|). 

Our finite-size estimates together with the extrapolated values for v L , z/" and x c are 
presented in table 1. One can notice the excellent agreement of these extrapolated values 
with the exact results v L = 1/2, x c = 1/3, while the present result v" = 0.81731(5) 
corroborates the best earlier estimates of the longitudinal size exponent. 
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4. 6*-Point behaviour 



4-1. Determination of the 9 critical temperature 



It was shown in |y| that for non-interacting directed animals the dependence of the 
average site density (|32]) on the strip width is of the form 



Pn (w = I) ~ C n—~ l (40) 

for large n (here and in the following we use C without super- or underscript as generic 
notation for an undetermined constant). By analogy with this behaviour and with 
the case of isotropic lattice animals p9f this density is expected to have the following 
finite-size scaling form in the vicinity of the IDA #-point: 

p n (w) ~ r^~ l H[{w - w{6)) n 1/u t], (41) 

for n ^> 1 and w — w(9) <C 1, with 

1> = {l-4)/»o, (42) 
where v\ and Ug denote respectively the longitudinal and perpendicular size critical 
exponents at the (9-point. H(z) is a non-singular scaling function and vjr, the 
perpendicular thermal correlation length exponent, is expected to be the same as the 
perpendicular exponent vjj P for directed percolation. Using ( |32| ) to calculate p n we have 
determined the quantities 

, / v log(p n /Pn+l) 1 

Wn{w) = - — + 1. 43 

log [n (n + 1) 



Due to the scaling form (|4l| ) one expects that ijj n (w(9)) — > ip for n ^> 1, so the solutions 
of the fixed-point equations 

ip n+1 (w) = ip n (w), (44) 

yield a sequence of estimates w n {6) and ^e, n for the critical weight w{6) and the ratio ip. 

The results are given in table 2. Note that this two-parameter phenomeno logical 
renormalization (PR) group method gives a rather accurate estimate for the Boltzmann 
critical weight, w(9) = 3.814 53(1). This estimate is even better than the one obtained 
using relation QT5D and the result of a one-parameter PR approach for the directed 
percolation threshold It is, however, less accurate than the best one: w(9) = 

3.814 5244, deduced from ( |T5] ) and the value p c = 0.644 700 185(5) recently determined 



using a one-parameter series expansion method for DP on the square lattice ||36|| . So, 
we shall use this last value of w(9) in order to estimate z/jj and . 

4-2. Geometrical critical exponents v\ and 

There are different methods to calculate finite-size estimates of two independent critical 
exponents v\ and . The simplest one is to suppose that the finite-size scaling form (|33|) 
holds for each fixed w, which allows one to apply the approach described in the preceding 
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Table 2. Finite-size estimates of the critical interaction w(6) and of the geometric 
exponents at the collapse transition. The estimates w n (0) and ipe,n are obtained 

for w = w(9) = 



from ( 44 ) , while Vg n and lA _ are computed using equations 



3.814 524 4. Extrapolation 1 is based on the sequential fit (|39|) restricted to two terms. 



n 


w n (U) 




,,-L 
v 0,n 


v 0,n 


4,5,6 


3.84492999 


0.765241160 


0.434318716 


0.673635137 


5,5,7 




(\ 7£fl1 1 770C: 

U. /oUll / iZo 


U.4ozozDDlt> 


U.D 1 004yzL)l 


6,7,8 


3.82611670 


0.756970160 


0.431222068 


0.676643208 


7,8,9 


3.82242111 


0.754906943 


0.430561122 


0.677317531 


8,9,10 


3.82001545 


0.753483968 


0.430142892 


0.677757109 


9,10,11 


3.81868653 


0.752461914 


0.429866869 


0.678056240 


10,11,12 


3.81769328 


0.751703212 


0.429678725 


0.678266797 


11,12,13 


3.81699662 


0.751124448 


0.429547246 


0.678419088 


12,13,14 


3.81649309 


0.750672740 


0.429453585 


0.678531701 


13,14,15 


3.81611976 


0.750313284 


0.429385894 


0.678616506 


14,15,16 


3.81583685 


0.750022425 


0.429336467 


0.678681344 


15,16,17 






0.429301202 


0.678731277 


Extrapolation 1 


3.81453(1) 


0.7479(1) 


0.4293(1) 


0.6789(1) 



subsection. Thus, using three-strip PR for w = w(6) we have obtained two sequences 
of estimates v$ n and v^ n . The results are shown in table |2|. 

A more sophisticated analysis consists in taking into account explicitly in the 
scaling function the two independent relevant fields that control the tricritical behaviour 
near the #-point. One then writes 

£ll ~ n K<, F[n 1/ui u, n l/u ' 2 v], (45) 
where to leading order the scaling fields are related to the model parameters by 

u = a(x — xg) + b{w — w{6)) + ■ ■ • , (46) 

v = c(x — xg) + d(w — w(6)) H , (47) 

and 

n e = vl/vfr. (48) 

The smaller scaling exponent, Vi, say, is related to the transverse size of the animals 
at the #-point PSI, so V\ = v^, while the larger one controls the thermal correlation 
length, so v 2 = . 

Now fixing w = w{6) one can still use ( P5j ) and to obtain a sequence of estimates 
xe, n and Kg in for xg and Kg respectively The extrapolated values for xg (see table ||) are 
in excellent agreement with the value deduced from ([14]), xg = 0.060 049 912 4(10). The 
main uncertainty comes in fact from the value of p c used to estimate w(6). 

Taking the partial derivatives with respect to x and w for a given width n (we 
follow closely the method and notations of p9|), one has at the ^-point 

a^jl = n K0 (a n 1/ui dF/du + c n 1/u2 dF/dv), (49) 
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ot^ll = n *e(b n 1/vi dF/du + d n 1/v2 dF/dv). (50) 

The unknown constants a dF/du and h dF/du may be eliminated using another width 
to: 

8 <) 



c (n 1/v2 - 1/vi - m 1/v2 ' 1/vi ) dF/dv, (51) 
d (n 1/va - 1/vi - m llv2 ' l/ui ) dF/dv. (52) 



The ratio c/d may in turn be eliminated using a third width I, yielding finally an 
equation for the exponent z = k$ + l/v\ : 

~ (n/m) z d x j\ = dJil - (n/m) z d w £j ^ 
d x £,l - (n/l) z d x Q d w & - (n/l) z d m Q 
The same calculation may be performed again exchanging the terms involving u% and 
z/ 2 . One finds that z' = hq + l/v 2 also satisfies (^3|), so one expects this equation 
to have two positive solutions, from which estimates for V\ and v 2 may be obtained 
since Kg has been estimated independently. The numerical results obtained using three 
consecutive widths are presented in table where it is readily seen that the extrapolated 
value Vi = 1.097(1) agrees with the value of the perpendicular percolation exponent 
v dp = 1-096 854(4) |3(|. The ratio Kg is also very close to its directed percolation 
counterpart z^p/^dp = L580 745(10) 

This suggests that the correspondence with critical directed percolation clusters 
is also valid for the geometric properties of IDA clusters in addition to their thermal 
critical behaviour. Indeed, one finds that ul/ujj P = 0.391 55(7) agrees with the value of 
the exponent a = 1/(Pdp + Jdp) = 0.391510(2) [06| . This is indirect analogy with the 
relation v c = a v v v for isotropic percolation, where u c is the exponent associated with 
the typical radius of large clusters p3| . Although interacting animals are related to 
directed bond percolation through ( |10"D it was not abvious a priori that their geometrical 
exponents should also be related simply to those of DP, as they depend on the behaviour 
of Gqr(x, w) at fixed w and not along Dhar's line. 

4-3. Density exponents 

Defining now the thermal density exponent f3 t through the variation of the density along 
the critical line : 

p{T) ~ C {w - w{0)Y\ (54) 
one obtains from the scaling relation (|4l|) and the numerical results above for ip and 

A =(1-^)^ = 0.2765(4), (55) 
in agreement with the percolation exponent (3 dp = 0.276 486(8) p6 |. 



At a fixed temperature, on the other hand, the density in the gel phase (for 
x > x c {w)) is expected to vary close to the critical line as 

p(x) ~C (x-Xciw)) 13 *, (56) 
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Table 3. Finite-size estimates for w — w(0) = 3.814 524 4. The estimates xg, n of the 
critical fugacity are obtained using the PR equation (|35|), while those of the critical 
exponents rely on ( |36"1) and ((5^). Extrapolation 1 is based on the sequential fit ( |39|) 
using only two terms (i.e., three parameters), except for the case of xg >n , where a 
five-parameter fit was used. As before, the BST algorithm is used in extrapolation 
procedure 2. 



n 


Xe,n 


t,II /»>-L 

v e,J u e,n 


V 0,n 


,,-L 
V t.n 


6,7,8 


0.0600496104928 


1.56917012 


0.431744944 


1.12856466 


7,8,9 


0.0600497568952 


1.57311880 


0.430976762 


1.12439401 


8,9,10 


0.0600498269551 


1.57570908 


0.430445690 


1.12138438 


9,10,11 


0.0600498630888 


1.57727568 


0.430203390 


1.11854334 


10,11,12 


0.0600498828566 


1.57854083 


0.429910536 


1.11679217 


11,12,13 


0.0600498941942 


1.57937827 


0.429743411 


1.11515184 


12,13,14 


0.0600499009511 


1.57998468 


0.429621850 


1.11377199 


13,14,15 


0.0600499051063 


1.58043100 


0.429531918 


1.11259417 


14,15,16 


0.0600499077280 


1.58076371 


0.429464479 


1.11157632 


15,16,17 


0.0600499094171 


1.58101415 


0.429413361 


1.11068746 


Extrapolation 1 


0.0600499124(2) 


1.5805(5) 


0.4294(1) 


1.097(1) 


Extrapolation 2 


0.0600499123(2) 


1.581(1) 


0.4293(2) 


1.097(1) 



where the exponent (3 X depends on the temperature and is related through scaling 
relations to the other critical exponents, yielding 

p x = u j + - 1 = 0.1082(1) for w = w e , (57) 

= J + p 1 - 1 = 0.3173(1) for w <w e . (58) 

We have not checked these predictions numerically. 

5. Low-temperature region 

5.1. First-order transition and essential singularity 

For w > w(6) typical animals are in a collapsed state and cover a finite fraction of the 
lattice sites in the thermodynamic limit. One expects the transition at fixed temperature 
to be first-order P9|| and the average site density p(x,w) to jump from zero to a finite 



value on the critical line given by x c = (w — 2)/w(w — l) 2 in that region. 
Using relation fllOD the animal generating function may be written 

in — 9 

G(x,w) = Z N (w)x N = (- ) J2 Pr ( N )- (59) 

TV W N 



If we now assume, as discussed in section ( p.3|) , that for all p > p c the probability of 
large finite directed percolation clusters varies as 

log Pr (N) ~ - C N c , (60) 
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Table 4. Finite-size estimates of the critical fugacity x n (w — 8), denned by 
Knax(xn) = lj and parallel correlation length on Dhar's line = — 1/log X n (xi}), 
with x D {w = 8) = 3/196. 



n 






x D 






4 


0.015306166798839312837 


4.43 


10" 


8 


8.8814536150 10 4 


5 


0.015306125866922874256 


3.42 


10" 


9 


9.2199061032 10 5 


6 


0.015306122730364804891 


2.81 


10" 


10 


9.3327985279 10 6 


7 


0.015306122473218244718 


2.42 


io- 


11 


9.2866324597 10 7 


8 


0.015306122451137529752 


2.16 


io- 


12 


9.1271646904 10 8 


9 


0.015306122449176587816 


1.97 


io- 


13 


8.8871999634 10 9 


10 


0.015306122448997933046 


1.83 


io- 


14 


8.5908642935 10 10 


11 


0.015306122448981326805 


1.73 


io- 


15 


8.2562064446 10 11 


12 


0.015306122448979758112 


1.66 


io- 


16 


7.8968815836 10 12 


13 


0.015306122448979607947 


1.61 


io- 


17 


7.5232899054 10 13 


14 


0.015306122448979593412 


1.58 


io- 


18 


7.1433739792 10 14 


DC 


0.015306122448979591837... 












where ( < 1 , equation ( |5T?D is similar to the one obtained by Fisher jE| for the droplet 
model of condensation, with ( corresponding to a surface energy exponent. He showed 
that G(x, w) can be expressed as 

G(x,w)= / f(t)dt, (61) 

Jo x c e l — x 

where the auxiliary function f(t) vanishes at the origin with an essential singularity of 
the form f(t) ~ exp (— C jt^ ), so G(x, w) as well as all its derivatives remain finite when 
x — > x c from below and G has a very weak singularity of the form 

C 

log G sing — — , (62) 

yX c X J 

with (' = (/(l — ()■ I n particular ^' = 1 if £ has the same value 1/2 as for isotropic 2d 
percolation. 

5.2. Numerical results 

As discussed in [B9], WM the phenomenological renormalization approach may still be 



used to determine the critical fugacity, as the parallel correlation length diverges at the 



transition. A direct application of the 3-strip PR approach using (^ - 138|) yields results 
for x c that converge very rapidly to the value expected from (|i~4l) . As can be seen from 
table § (see the third column), the finite-size estimates x n of x c are already quite good 



for narrow strips. To obtain such an accurate numerical check of (14) we had to use a 
mathematical package allowing for high-precision numbers (the precision of our input 
data was set up to 30 digits). 

The exponent K n , however, is found to diverge, signalling a failure of ([33]) to 
describe adequately the finite-size scaling of A more detailed numerical analysis 
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of the behaviour of for an interaction w = 8 deep in the low-T region, shows that it 
increases exponentially with the width n if one fixes x at the critical value x c (w) given 
by (0): 

£l(x B ) ~ AeMKn), (63) 

and that the fugacity x n at which £| diverges (i.e., A n (5: ri ) = 1) converges exponentially 
towards x c (see table These results suggest the following finite-size form: 

C C 

/HI ~ ~ (Q4) 

W (x„-x) A (x c -x + D exp(-K'n)) x ' 1 ; 

where A = K/K'. Numerically one finds that K ~ K' = 2.21(2), so we conjecture that 
A = 1. 

The average mass of the clusters 

< N(x) > = x — ^— (65) 
ox 

is dominated by the regular part of G(x) and remains finite when x — > x c from below at 
fixed temperature, as well as all the higher moments of the mass distribution, so the mass 
and linear dimensions of the typical clusters also remain finite. This implies that the 
characteristic transverse length to be compared to the strip width in a scaling expression 
remains finite at the transition, which is consistent with the observed divergence of the 
effective exponent K n . 

In the zero-temperature limit (w — > oo) the critical fugacity x c — > 0, with x c w 2 — > 1, 
and the generating function G(x, w) is dominated by the configurations with the largest 
possible number of pairs of neighbouring occupied sites at fixed N. These correspond 
to compact animals, having no internal holes and the maximum number of independent 
loops, or equivalently the minimum number of lateral surface sites (i.e., sites with either 
the smallest abscissa on a given row or the smallest ordinate on a given column of the 
underlying square lattice). One would expect these animals to be globular rather than 
elongated, but it is not clear to us how to extract from the transfer matrix results values 
for the geometrical exponents of the (extremely rare) large clusters in that regime. 

Let us note that the present compact animals are different from the "directed 
compact percolation clusters" studied by Essam HS[ , which just obey the constraint that 



two diagonally neighbouring sites are necessarily followed by an occupied site. These 
are themselves a subset of another class of compact clusters studied in the literature, 
named "fully directed compact" by physicists |R| and more precisely "diagonally 



convex directed" (DCD) by mathematicians [f48~j |49| . DCD animals may be characterized 
as having the minimum accessible perimeter for a given length (i.e., unoccupied sites 
having an occupied predecessor). They have been exactly enumerated with respect to 
various parameters HH, |50|, f[l[] and are expected to be extremely elongated, with v" = 1 
and v L = (where i >± is associated with the local width at a given large distance from 
the origin, rather than to the overall animal width 
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0.40-n 1 1 1 1 1 1 1 1 1 1 r 




Figure 2. Critical line x c w versus interaction parameter w. The part of the solid 
curve for w > w{9) corresponds to a hrst-order transition and is given by (|l4|), the 
dash-dotted line is its analytic continuation for w < w{9). The symbols denote special 
values of w at which G(x, w) is related to the generating functions of various simple 
ensembles. 



6. Special points 

As pointed out in section (|2.4j ), the phase diagram of interacting directed animals in 
the (w,x) plane contains several special points related to various ensembles of animals 
having specific geometric properties. 

• Thus, from ( p2]) the generating function of bond animals is obtained on the line 
w = (i-x) • Using strips of width (14, 15, 16) we find this line intersects the 
critical line x c (w) for w c = 2.285 08(1). This is in agreement with the value 
yb = 0.285 0875(8), obtained for their critical fugacity by Conway et al. since 
in the present notations w c = 2 + yb. 

• Bond trees correspond to w = 2 according to (|24D and we find that the 
corresponding x c = 0.147892 24(1), in agreement with (but more precise than) 
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the critical fugacity yu = 2x c = 0.295 785(10) found in pS . 



Site trees correspond to the limit w — > with y = fixed and finite. They 



were studied as Model B in |15| and their critical fugacity was found there to be 
y st = 0.368 649(2). 

• Finally for a site fugacity x = 1 one recovers the generating function of strongly 
embeddable animals, see (^7j ), so their critical fugacity y em b = w c (l). We find 
Vemb — 0.356 563 26(1) for this quantity, which does not seem to have been estimated 
previously in the literature. 

A phase diagram summarizing these various results is displayed in figure |2|. It is 
noteworthy that for w < wg the plotted quantity y c = wx c (w), which according to ([18|) 
is just the critical fugacity for animals enumerated according to loop number, varies 
quasi linearly with w, while this is not at all the case for the analytic continuation of 
the low-temperature expression (|i"4|). 



7. Crossover region close to the collapse transition 

7.1. Critical line and crossover exponent 

We now focus our attention on the region close to the #-point. As noted above, for 
w > w(6) the critical line x c (w) is given by the simple expression (|i~4l) and a remarkable 
consequence |L8| is the absence of a singularity in the bulk free energy per site f(T) 
when w — > w(8) + . The collapse transition is just signalled by a singularity in the size- 
dependent correction term, which may be interpreted as the vanishing of a surface free 
energy. For w < w(8) the critical line departs from the analytic continuation of (|14|), 
which we denote x^ in the following to avoid confusion. The singular part is traditionally 
written under the form 

x c - x D ~ C(w(9) - w) 1 ^, (66) 

where is the crossover exponent. 

In most cases it is quite difficult to obtain accurate values for such crossover 
exponents from numerical data, as generally neither the non-singular part of the critical 
line nor the value of the critical temperature are known with high precision, making a 
fit of the form flBTf ) subject to large uncertainties. Here the situation is quite favourable 
as we can take advantage of the analytic nature of the free energy in the low-T phase 
and of the high-precision knowledge of the critical point w(9). The remaining difficulty 
comes from the uncertainties due to the extrapolation procedure used to obtain x c (w), 
which prevent us from analyzing reliably the data extremely close to w{6). Fortunately 
there exists an interval of values of w where the extrapolation uncertainties and the 
corrections to (|66"D are both small enough to allow a reliable analysis. 

The numerical results are displayed as a log-log plot in figure |3| and the slope of the 
curve yields 

1/0 = 2.53 ±0.05 (67) 
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log[w(e)-w] 



Figure 3. Singular part of the critical fugacity, versus distance to the collapse 
transition, in a doubly logarithmic plot. The line is a linear fit to the numerical 
results (filled squares), whose slope yields the value = 2.53(5) for the crossover 

exponent. 



(j) = 0.395 ± 0.008. (68) 

This corresponds to a weak singularity for the specific heat measured along the critical 
line, as the associated exponent a |k| is negative: 

a = 2- 1/0= -0.53 ±0.05, (69) 

in strong contrast with what is found for the #-point of isotropic animals, where 
a ~ 0.48 |2"5fl. This also implies that in the vicinity of their collapse transition interacting 



directed animals should display geometric properties close to those of critical directed 
percolation clusters up to large scales. 

7.2. Relation with a percolation exponent 

In order to derive a relation between the crossover exponent and a directed percolation 
exponent, we rewrite ( |T2"D under the form 

logx c = logx D - lim (1/JV) log Pr(iV), (70) 

N— >oo 
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where xjj = (w — 2)/w(w — l) 2 . For w < w(9), the bond probability p < p c and the 
cluster probability Pr(N) is expected to decay exponentially with N, so the second 
term on the right-hand side of (70) is no longer negligible when iV — > oo. Equation (|70|) 



may then be written 

(1/N) log Pr (N) ~ -log(x c /x D ) for N > 1, (71) 
~ _ _5 £ f or x c — x D <tix D , (72) 

~ C[w(0) -w] 1 /*, (73) 

using ( p^| ) in the last line. Relation (|TT|) between p and w still holds in the high-T phase, 
so w(9) — w ~ (p c — p)/(l — p 2 ) an d we finally obtain 

(1/N) \ogPr(N) ~ C(p c - p) 1 ^. (74) 

Now, close to p c the probability distribution of cluster sizes is expected to follow a 
scaling law of the form |54] 



Pr (iV) ~ N 1 ^ ^[N a (p c - p)], (75) 

for N ^> 1 and 2 = N a (p c — p) finite, where ^(z) is a scaling function finite for z = 
and we have followed the notations of |23j for the Fisher exponents r and a. 
For these two expressions for Pr(N) to be compatible it is necessary that 

<p = a. (76) 

A very accurate value for a may be deduced from Jensen's recent results [|36| and a 
standard scaling relation between percolation exponents : 

* = 1/{Pdp + 1dp) = 0.391 510(2), (77) 

in good agreement with our numerical value (|68D for <p, which was obtained without 
invoking scaling assumptions. 

8. Conclusion 

We have shown that the combined use of the transfer matrix method and of a 
phenomenological renormalization-group analysis gives very accurate results when 
applied to interacting directed animals (IDA), in particular close to their collapse 
transition. The present results are still not quite as precise for the tricritical fugacity and 
the tricritical exponents as those obtained from the best series calculations for directed 



percolation [|6|, but the method is more flexible so it could also be applied without 
excessive effort to situations involving more complicated interactions. 
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